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as INTRODUCTION 


The ccmplexity involved in designing a control system 
cf reasonably low crder for dynamic systems cf substantially 
higher order is a subject about which much has been written. 
The list cf references at the end of this study include only 
a small ferticn of the material available concerning this 
subject. A néew development has been intrcduced by Maurice 


F. Hutton and Bernard Friedland, [Ref 1], which provides a 


systematic and analytical approach =O obtaining 
reduced-crder approximations from higher-crder transfer 
functions. This method is referred to as the "Routh 


Approximatice Methcd," and is based on an expansion that 
uses the Rceuth table of the original higher-order transfer 
function. Tke desire for low ocrder models to  Sirulate 
higher-order practical systems, such as electrical fower 
plants, chemrical processes, aircraft model designs and 
electrcnic cirfcvitry is obvious when dealing with such 
ccmplex systezs. The "Routh Approximation Method", is thus, 
an extremely powerful tool to utilize in dealing with these 


types of proklems. 


In ccmpariscn to other available methods fresently 
used, the Kotth methcd has many advantages. The widely used 
"“NDcminant Pcle'', approximation method , which is based on 
approximating a system by ete el ez the poles, Or 
characteristic rocts nearest the imaginary axis, has the 
disadvantage that the roots of the characteristic equation 
must be found. For a very~high-order system, this is not a 
trivial task. The "Pad‘'e" approximation is tased on setting 
the numerator and denominator crders tc a desired value, the 
ccefficierts are then chosen so that the Taylor series 


expansions cf the approximant and the original transfer 
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function agreé€é in aS Many terms aS possible. This method 
Ereduces accurate results, however, it is limited in 
applicaticn tc single-input, single-output system analysis. 
In additicn, an unrstable approximant may ce obtained frem a 
stacle system since the approximant's poles depend on toth 


tke original eguation's numerator and denominator. 


The "Reuth Approximation Method," preserves stability 
if the criginal transfer function is stable. It provides an 
efficient reans of obtaining lower-order approximants for 
Multiple inputs or ottputs and is very adaptable to computer 


programming. 


This thesis 1s concerned with the development of a 
computer frcegram which utilizes the "Routh Approxigation 
Method," to obtain lower-crder transfer functicns from 
higher-crder functions and compares the resulting éguations 
to the cridgiral eguations by graphically displaying the 
response cf each tc various input excitations. The prcegran 
was designed to provide output data which is useful in 
determining which degree of approximant is best suited to 
Simulate tke original higher-order equation. The crders 
available are the first through the fourth order reduced 
transfer functicns. The program is capable of rec¢ucing 
transfer functions up to, and including the tenth orcer and, 
Without loss of generality, may be extended to handle an 


almest unlimited order. 


For iliustrative purposes, a seventh crder system was 
used, aS an exampie, to demonstrate the computer frogranm's 
Capabilities. This particular system was chosen merely to 
demcenstrate the simplicity involved in utilizing the prcegranm 


and has ne physical relationship to any real systen. 


il: 


eee ieocunee cr THE PROBL EM 


In @ppreximating higher-crder systems by lower-order 
models, a linear system is desirable since the linear 
characteristic equations are less complicated than the 
non-linear eégquaticns. Therefore, this study 1s primarily 
ccncerned with linear, time invariant characteristic 
equations, since the objective in develcfing the conputer 
prcgram sas te illustrate a simple analytical approach to 
obtaining lcwer crder characteristic eguations frem higher 


order systers. 


oc is obvicus that a reduced-order model cannot 
characterize a given system as accurately as a higher-order 
model. The validity of the lower-order model is based upon 
its degree of succeSS in approximating the higher order 
system in representing the characteristics of primary 


Mmircerest. 


Interpreting the soluticn of a higher-order systen 
cften results in computational difficulties which are 
recuced Ey appropriate selection of a reduced-order 


apPEICcXximaticn. 


Ideally, a reduced-order model would approximate the 
higher-orcer system in both low and high frequency ranges. 
In deing sc, Scme accuracy is lost in compensating fcr the 
different responses of the systen wile Variations in 


frequency. 
The lew frequency model more closely approximates the 


higher-order system than a model composed cf both lew and 


high crder frequency characteristics. The procedure 


1 


utilized in this study places emphasis on the low freguency 


model. 


For unstable systems, the program described herin is 
still valid, however, the original higher-order systen must 
first be mcdified by shifting the imaginary axis prior to 
the apprceximation. Thus (EO an unstable transfer 
function,#(s), the equation must be changed tc H(sta), where 
a is chcser sufficiently large so that H (s) is 
asymptotically stable. This procedure is described in 
G@etail itr fkef t], pp 332. 


LS 


MGieecoNeUrha PROGRAM CRITERIA 


A. GHNERAL 


The frogram, referred to as ROUTH1, waS written with 


meme £eclicwing criteria: 


ROUTH1 consists of less than 150K of storage and taxes 
less than 2 minutes of computer time. Tiers. nC ca One? 
conserves efficiency, but also provides the user with the 


desired data at a minimum cost. 


2. Ease ci use 


ee 


Input data reguired consists of eleven data cards for 
Maximum utilizaticn of the program. Emphasis is placed on 
the ¢€ase cf using the program to obtain the desired results 


With minimum time expended on computer programing. 


a. DiometacCst second, third and fourth order 
approximants to the original equation are printed in 
transfer function format. Both numerator and denoninator 


coefficients are printed in ascending powers of S. 
pee The roots of the original and reduced equations 


are previded to enable the user to study the response cf the 


systems in the frequency domain. 
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¢. Choice of variables for print out in table forn 
is available witn up to eight variables maximum for any one 


cf three fessible runs. 


diz Choice of variables for graphical output with up 
to four curves may be plotted separately cr all cn one 


gaa Oh . 


CeeeGEdvaecdl  OUEpUE Eesponse tS input excitations 
corsisting Ciena mnotece name, ©L Sinusoidal input are 
available. fhe graphs display the original output response 
compared to the lower order response and the error is 
pictted tc display the differences in response. Muitiple 


inputs may be used. 


Dee SUBRCUTINES 


Imm StrrOutines are utilized in ROUTH1. The rocts of 
the original and reduced equations are determined by 
Subroutine PRQD , waoaich was taken from [Ref 3], and the 
tables and graphs are determined by subroutine RSDUCT1 which 
fomeae modi tication Of INTEG? ; {Ref 3]. 
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IV» FOURTH ORDER EXAMPLE 


Pema oft P PROCEDURE 


To illustrate the "Routh Approximation Method", A 


fourth order example was chosen for simplicity. 


The method, described in detail in [Ref 1], consists of 


three basic steps. The following transfer function is used 


to illustrate the procedure: 


H({s) = 1 


20 + 328 + 24% 4+ 384 ¢ 


The first step is to compute what are termed Alpha and 


Beta coefficients from the Routh Table shown on page 19. 


— 


Alpha, =0.625 is; I 
Alphaj=1.684 6.3158 
Alpha.,=3.008 ib 
Alpha ,=6.3158 

Beta i 


Beta, =0.03125 = 05.25 
Beta,=0.0 0.0 
Beta.=0.03958 

Beta ,=0.0 
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Step two in the procedure is to obtain what are termed 
the Routh convergents, which are based on the following: 
Letting A, (s) and B, (s) denote the denominator and the 


numerator, respectively, of the ee Routh convergent, 1.ée., 


A,(s) = Alpha,s + 1 

B, (s) = Beta, 

A, (Ss) = Alpha, Alpha, s* + Alphas se iL 
B, (s) = Alpha,Beta,S + Beta, 


The general expression from [Ref 1] is the following: 


A, (Ss) = Alpha, (s)A,_, (s) + A, 9 (S) 
B,. (s) = Alpha, SB, _, (Ss) + By. (Ss) + Beta, eae, 3 sans 
with A_, (s)=0.0 B_, (s)=0.0 
Ay(s) =1.0 By (s) =0.0 


The Routh convergents for the fourth order example are 


the following: 
R, (s)= Oy 6) bras 
O Meese os Al 


R,(s)= O05 6526 


1.05263S* + 1.68425 + 1.0 


a, ee 0.158833S* — 0.0083 


Goan MalGea sn Ge6cas + 1 


R,(s)= s? 


2087 - 32S° oo 249° + 8S + 1.0 


ny 


The third, and final step in the procedure is to apply 


what is termed a reciprocal transformation, defined by 


H,(s) = =xR, (1/s) 


which 1s merely a reversal of the order of the polynomial 
coefficients. Thus, for the example given, the reduced order 
approximations are given by the following: 


H, (s) = OO S125 


a oS 
H,(s) = 0.056316 

s* + 1.68425 + 1.05263 
H,(s) = -0.0083S* + 0.158333 

s? + 3.633S8° + 5.0665 + 3.166 
H,(s) = i 


s* i gs? +2487 aoe a 20 


As expected, the fourth-order approximation is the same 
as the original equation. 

The poles of the approximants are illustrated in table 
IV.2. The poles of the approximants approach the dominant 
poles of the original higher-order equation as the order of 


the approximant is increased. 
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Veeecourlurar SAAMPLE 


Phe oeWee NTH CRDER SYSTEH 


A seventh order system with the transfer function given 


by [Ref 2], as 


384x10/ 


G(s)=— Se ae a 5 ao 
s/44325°4626708>+3615900S'+75114000$°+5539200008-... 


..+1443200000S+384x10/ 


Spel CLaccOoread tom as 


384x10/ 


G(s)= 
($*4+2S+10) (S+10) (S+20) (S+80) (S+120) (S+200) 


was reduced to the low order approximants Dy 
utilizing ROUTH1 and the response to input excitations, 


consisting of a Step and Ramp are illustrated in figures 5.1 


wmanmouga 5.6. 


The rocts of the system and its low order approximants 


pec atliuetrated in table V¥.1. 
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Poles and Zeros of Approximant 


-1.0 + 32.999 
R, |-10.89 + 32.30 {+ 316.58 


2 
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-2 
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POLES OF ROUTH APPROXIMANTS 
Table V.L 
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VeweeecoreannosON LTO OTHER METACDS 


A. DOMINANT FOLE MeTHOD 


The Deninant Pole Approximation method 1S paseG upon 
utilizing the poles closest to the imaginary axis. Tne 
eguation must be factored to obtain the characteristic 
rocts. Fer the given seventh-order system, the lower- order 


equations are given by this method as the following: 


10 

CS) aoe 

S$~+2S+10 

smeye 

1) 

S~+14S5°~+30S5+100 

2000 

SS) ar a 

S°+34S° +310S +700S+2000 

Graphical flots of the Dominant Pole reduced 


equations, in response to Step inputs, are illustrated in 
muiceaeo,1  thEOUugh 6.3 in comparison to the original 
equation and fegquecs@o.st throtgh 6,6 illustrate tne 
ccmpariscn tc the Routh equations. The error between the 
Systems is also plotted. Table VI.1 gives the analytical 


data for perficrmance measure comparisons. 
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Peek At EVES OPLTEUETZATION METHOD 


The fcllowing equations, representing the seventh-order 
example, were taken from [Ref 2]. This method maxes use of 
an iterative finimization technique to locate the best fole 


and zero locations for the lower-order nodels. 


ea O5135 6 


SS a 
Stee colGed, «2055516 


8) an Pas aud 


S 2 


dS) 2-3 
meee oo oat hoe 0 lSSto2. 286) 


1470.1403 


~ 3 Z 


1 (s) = So 
Seemeerz04s +2082 58425 +552.5241541470.1403 


Plots of these equations versus the Original 
equation and the Routh equations are illustrated in figures 
6.7 through 6.12 in response to Step inputs. Table VI.1 


provides numerical data for performance measure compariscnas. 
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Seer nheorNTATION OF DATA 


Comparisons were nade between the "Routh Approximation 
Method", and that of the Dominant Pole and Iterative 


Optimization methods previously described. 


iiembaslsesOm COMDadrison consists of the following 


materia : 


1. Peak ONISEES ITN BS Ss WE 


Max a) 


ZS hay Liane -———_— = T,=time PommrGt) tO reach 0.5 Y 


the very first time. 


S 


Sere Soe ne] —-—-———— T =time Mom) —eOug@uncomn, O.1, Lo 
0.9 of the final value. TW=1/5W. 


4. Settling Time---—-T_=time at which Y(t) =¥_ 
5. Graphical Representation in response to Step inputs. 


6. Average Error----J= DE/t.,| ; t,=integration steps 


Table ¥VI.1 illustrates the above comparisons. 
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SYSTEM 
TYPE Mpt 


O-4 | Le 205 E 533 re fatal 5.866 
3 
| D-4 | Ry: I 528 ‘0 qa is 866 1 
| . 22 ee eee 
R-3 14.247 10.533 
: { 
} ee 
| ! : 
0-3 1.320 10.533 
| ne 
i 
p-3. | 1.346 [0.502 10.711 [5.866 1.15x10 at 
| 1 | = a foes m! : 
pe | 1.083 10.400 9) 7k 15.866 0.062 | 
| ) ae ee ee aa 
| | 
0-2 {1.285 40.533 Do pal 5.866 a. : 
| ; | : 
D-2 [1.350 )0-412 10.684 |5.866 10.059 : 
TUNEHES | WSEey al 
Table A = Original Seventh Order 
symbols R = Routh Aporoximation 
D = Dominant Pole Approximation 
O = Optimum Minimization Method 


PERFORMANCE MEASURE COMPARISONS 
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Vier oONe LUST ON 


The grafhicai data and tabulated analysis of this 
thesis irdicate that the Routh Approximaticn Methced is a 
valuatle formulation technique which produces very 
satisfactcry results in acquiring approximations to 


higher-order systems with minimum cost. 


In ccmpériscn to the other methods cf analysis, the 
Routh Arpreximation Method offers a quick and easy 
analytical approach to obtaining low order models. Unlike 
the Fad'te Arpreroximation, the Routh method ensures statility 
of the lcwer-crder models, if the original higher-crder 
system is statle. The original system need net be factored, 


as in the Dcninant Pcle Approximation method. 


The ccaputer program, ROUTH1, utilized to acquire the 
recuced crder equations, the roots cf the eguations and _ the 
graphical ficts and numerical tables, takes considerably 
less time than the minimization technique utilized in 
reference 2. However, the minimization téchnigue operates 
efficiently withecut prior knowledge of the system's transfer 
mane clon . In the Routh method the availability cf the 


transfer function is a necessary requirement. 


In ccmparing the low order equaticns of the various 
methcds discussed, the Routh method has proven to be a valid 
and efficient solution to the problem of cbtaining gocd low 


crder approxigzants to complex higher-order systems. 
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